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Abstract 

A procedure for designing a fuselage having a 
prescribed effective area distribution computed from 
—90° Mach slices is described. This type of calcu- 
lation is an essential tool in designing a complete 
configuration with an effective area distribution that 
corresponds to a desired sonic boom signature shape. 
Sample calculations are given for Mach 2 and Mach 3 
designs. The examples include fuselages constrained 
to have circular cross sections and fuselages having 
cross sections of arbitrary shape. For a prescribed 
effective area distribution having sharp variations, 
the iterative procedure converges to a smoothed ap- 
proximation to the prescribed distribution. For a 
smooth prescribed area distribution, the solution is 
not unique. 

Introduction 

One approach to minimizing the effect of sonic 
boom noise is to attempt to design the aircraft con- 
figuration so that its ground level signature has a 
“low boom” shape. A class of such low-shock and 
low-maximum-overpressure signatures has been stud- 
ied in references 1 and 2. Reference 2 also includes a 
systematic procedure for generating the effective area 
distribution for a configuration that produces the de- 
sired ground level signature along the flight ground 
track. This effective area distribution includes the 
distribution of areas of projected —90° Mach slice 
cuts (cuts by planes swept at the Mach angle and 
perpendicular to the x 2 -plane) through the config- 
uration. It also includes a significant contribution 
determined from the longitudinal lift distribution of 
the aircraft and a contribution due to nacelle-wing 
interference. 

In order to design a configuration that has the de- 
sired effective area distribution, one procedure that 
has been adopted is to start with an aerodynamically 
feasible configuration, approximate its total effective 
area distribution with a low-boom distribution, then 
redesign the fuselage only so that the desired distri- 
bution is obtained exactly. The problem is thereby 
reduced to that of designing a fuselage having a pre- 
scribed effective area distribution. This problem is 
treated in the present investigation in a general con- 
text. The fuselage may be cambered, it may be con- 
strained to have circular cross-section shapes, or it 
may be allowed to vary in cross-section shape. 

Symbols 

A(x) effective area distribution; area 

of projection on x == 0 plane of 
sector cut by —90° Mach slice 
that intersects z = 0 plane at x 


Ct ^-coordinate of fuselage camber 

line at terminal point 

lb body length 

l p length of prescribed area 

distribution 

M flight Mach number 

p iteration parameter 

r body radius for bodies having 

circular cross sections 

x, y , 2 Cartesian coordinates; x-axis 

taken through fuselage nose in 
flight direction 

P = y/M 2 - 1 

Subscripts: 

a average of maximum and mini- 

mum values 

d design 

i iteration number 

s slicing plane 

0, 1 0th and 1st iteration, respectively 

Procedure and Sample Calculations 

General Considerations 

For this investigation, the general approach is to 
gradually, by iteration, correct the fuselage shape to 
yield the desired effective area distribution, using a 
known analysis method. The fuselage equivalent area 
distribution calculations utilized a Langley version of 
a code originally written to compute wave drag from 
volume-effect Mach-sliced area distributions (ref. 3). 
The effective area associated with a particular Mach 
slice is assigned at the x-location where the slicing 
plane intersects the x-axis (see fig. 1). This effective, 
or projected, area is calculated by computing the 
coordinates of points of intersection of the slicing 
plane with the body shell, setting the x-coordinates 
equal to zero, and then integrating the resulting 
closed contour in the x = 0 plane. This projected 
area is times the actual area cut by the Mach 
slice. 

Circular Cross-Section Fuselage 

If the fuselage is cambered, cross sections are 
usually defined, for the sake of convenience, in planes 
perpendicular to the x-axis (x = Constant planes) 
rather than in planes perpendicular to the camber 
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line. For some designs these cross-section shapes are 
constrained to be circular. In this case, a change in 
the radius at any station affects the Mach-sliced areas 
cut by all slicing planes that intersect that circular 
section (fig. 2). Thus, it is not possible to change 
the effective area distribution of the fuselage locally, 
that is, at only one point. However, if the specified 
area distribution is smooth with gradual variations, 
it is generally possible to design the corresponding 
fuselage while maintaining the circular cross-section 
constraint. 

For such a fuselage, the design calculation pro- 
ceeds according to the following steps. 

1. Begin with an initial trial configuration having 
an effective area distribution that roughly approxi- 
mates the desired distribution. This is usually taken 
to be a body that has a length 

h = Ip + fiCt ( 1 ) 

where l p is the final value of x for which the effective 
area is prescribed, and Ct is the ^-coordinate of the 
camber line at the terminal point. The radii for this 
initial body are then computed by the formula 



2. Compute the actual effective area distribution 
Aq for this initial body, and the distribution of the 
error 

A,4(a;) = A d (x) - Aq(x) (3) 

3. For each Mach-slicing plane, determine the 
minimum and maximum x-locations of the intersec- 
tion of the Mach plane with the fuselage. 

4. For each input x-station within this interval, 
change the corresponding radius in such a way that 
the error A A is reduced. 

5. Iterate step 4 until the desired effective area 
distribution is attained. 

The crucial part of this procedure is step 4. As 
was mentioned earlier, the problem is that changing 
the radius at any input location changes the effective 
area for every Mach slice that intersects that circular 
cross section. For low supersonic Mach numbers — 
for which the slicing planes are nearly vertical — this 
problem is avoided if the input x-locations are suf- 
ficiently sparse. For higher Mach numbers, we rec- 
ognize the effect of overlapping intervals but rely on 
iteration and smoothing to yield the desired result. 

We denote the initial radius distribution in equa- 
tion (2) by ro(x). The new effective area distribution 
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Ai(x) should vary approximately according to 

r li x ) _ [A l( x ) , 4 x 

r 0 (x) Y A 0 (x) 

Now if Ai(x) is taken to be the desired distribu- 
tion A^(x), equation 4 should yield a radius distri- 
bution closer to the required distribution rg(x). An 
iteration equation is obtained from equation 4 by first 
subtracting 1 from both sides, 



then inserting an iteration parameter p, multiplying 
through by tq, and switching to iteration indices: 



(6) 


The iteration parameter p is used to control the 
size of the step taken at each iteration. It can be 
permitted to vary with x and/or with iteration num- 
ber, but the sample calculations in this study used 
constant values for p. A large value of p accelerates 
the iteration process, but too large a value may re- 
sult in an instability in the iteration. Consequently, 
small values (0.025-0.1) are often used, but for some 
relatively low Mach number cases, values as large as 
2.5 have been used effectively. 

Sample Calculations 

The computer code that implements this proce- 
dure utilized iteratively the analysis procedure of ref- 
erence 3. Its use is illustrated by the example shown 
in figure 3. Figure 3(a) shows the prescribed effective 
area distribution, as well as the effective area dis- 
tribution for the initial shape computed from equa- 
tion (2). Also shown are the results for the final 
design. 

In all the sample calculations, there is some devi- 
ation of the effective area distribution from the spec- 
ified distribution near the aft end of the fuselage. 
This deviation occurs because the value of A^i is 
not permitted to fall below a certain positive value 
since it occurs in the denominator in equation (6), 
and consequently, the iteration cannot converge in 
this region. 

Figure 3(b) shows the fuselage shape resulting 
from the design calculation of figure 3(a). Figure 3(c) 
shows the design with the vertical coordinates ex- 
panded by a factor of 6 to emphasize the variations. 

To illustrate the procedure for an M = 3 condi- 
tion, the same design distribution A d {x) was taken 
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as for the M = 2 case. Figure 4(a) compares the ob- 
tained effective areas with the prescribed distribution 
A d (x). The resulting shape is shown in figure 4(b). 
Figure 4(c) shows the design with the vertical coor- 
dinates stretched by a factor of 6. 

Smoothing, Solution Existence, and 
Uniqueness 

The computer code contains a provision for 
smoothing the fuselage after each iteration. However, 
little, if any, smoothing is actually required because 
the overlapping of intervals cut by different Mach 
slices results in an inherent smoothing. This over- 
lapping increases with Mach number, and so the con- 
verged area distribution will be quite smooth. Con- 
sequently, greater smoothness in A^x) is required 
at the higher Mach numbers. In fact, if Ad(x) con- 
tains corners or regions of very high curvature, the 
iteration cannot converge to A^x). In this case, the 
design calculation results in a smoothed approxima- 
tion to Ad(x). 

In an attempt to obtain effective area distribu- 
tions with regions of high curvature, a slightly differ- 
ent procedure was used. Instead of altering all the 
radii intercepted by a Mach slice, only the radii aft 
of the intersection of the Mach slice with the camber 
line were altered. This technique eliminated part of 
the overlap, and consequently the iteration converged 
rapidly except at very near the aft end. 

However, for the same input variables as for 
the previous example, the fuselage had a different 
shape — one with cyclic radii variations (fig. 5). This 
waviness occurs because the modified downstream 
fraction of the interval cut by each Mach slice, in 
attempting to correct for the entire interval, exag- 
gerates the correction locally. For example, if a 
particular Mach-sliced area is considerably less than 
that specified, the radii for the downstream part of 
the Mach-sliced interval would be increased signifi- 
cantly, because no contribution to the increase comes 
from the upstream radii. Now, however, downstream 
Mach slices that intersect this greatly enlarged part 
of the body have too great an effective area. The 
required decrease is accomplished by decreasing the 
downstream radii in an exaggerated manner, and 
so this cyclic behavior is propagated downstream. 
However, the desired effective area distribution is 
obtained. 

Thus, it is seen that, at least for some cases, 
the fuselage corresponding to a given effective area 
distribution is not unique. 

Arbitrary Cross-Section Fuselages 

When the cross section of the fuselage is not con- 
strained to be circular, a similar iterative procedure 


is used, but the design procedure is more straight- 
forward. In theory, one needs only to adjust the co- 
ordinates along each Mach slice in order to correct 
the effective area for that slice. However, the dis- 
crete input coordinates are located on x — Constant 
stations and not along Mach slice sections. Inter- 
polation could be used to compute the changes at 
each of these coordinate points. However, since the 
calculation must be iterated in any case, a simpler 
method was used. The intersection of the Mach slice 
with each lofting line was located, and the coordi- 
nates were adjusted at the input point on the lofting 
line that was nearest this intersection. (The jth loft- 
ing line connects the jth points on each of the cross 
sections.) 

The z-coordinates are adjusted relative to the 
central value of z 

^min («e) 


These values of z a (x) represent, to a close approxi- 
mation, the fuselage camber line. The iteration equa- 
tion (6) is now replaced by two equations: 


Zi{x) 


z a {x) + (Zi(x) - z a {x)) 


1 + 


Vi(x) = Vi - lW 


1 +p 






(8a) 

(8b) 


This iteration converges more rapidly than the circu- 
lar cross-section calculation. A sample calculation is 
shown in figure 6 for an M = 2 body. 

One problem that can arise with the use of equa- 
tion (8a) is that the maximum and minimum z- 
coordinates at each x-station are gradually changed 
in the design process. Consequently, their average z a 
is changed, and the redesigned fuselage has a camber 
line that is slightly different from that of the origi- 
nal fuselage. This variation is normally slight and in 
most cases does not represent a problem. 

However, the problem is avoided altogether if 
the z-coordinates are held constant and only the y- 
coordinates are iterated (eq. (8b)). Figure 7 shows 
the results of a calculation carried out in this way. 
Comparing this fuselage with that of figure 6 reveals 
somewhat greater lateral variation of the lofting lines. 

Concluding Remarks 

A procedure for designing a fuselage having a 
prescribed effective area distribution computed from 
-90° Mach slices has been described. This type of 
calculation is an essential tool in designing a com- 
plete configuration with an effective area distribution 
that corresponds to a desired sonic boom signature 
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shape. Sample calculations were given for Mach 2 
and Mach 3 designs. Examples included fuselages 
constrained to have circular cross sections and fuse- 
lages having cross sections of arbitrary shape. For 
a prescribed effective area distribution having sharp 
variations, the iterative procedure converges to a 
smoothed approximation to the prescribed distribu- 
tion. For a smooth prescribed area distribution, the 
solution is not unique. That is, more than one fuse- 
lage shape may have the prescribed effective area dis- 
tribution. This lack of uniqueness may permit some 
freedom in allowing for other design constraints, such 
as minimizing wave drag. 


NASA Langley Research Center 

Hampton, VA 23665-5225 

January 9, 1990 
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(a) Effective area distributions. 



(b) Designed circular cross-section fuselage. 
Figure 4. Example calculations; M — 3 design. 
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(a) Actual scale. 


Top 




Side 

(b) Expanded y- and z-scales. 

Figure 7. Fuselage designed for same conditions as for figure 6 but with j/- variation only. 
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